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Abstract 

The  forces  and  moments  produced  by  a  vehicle’s  aerodynamic  control  surfaces  are  often  nonlinear 
functions  of  control  surface  deflection.  This  phenomenon  limits  the  accuracy  of  state-of-the-art  control 
allocation  algorithms  since  all  of  the  approaches  are  based  on  the  assumption  that  the  control  variable  rates 
are  linear  functions  of  the  surface  deflections  and  that  control  variable  rate  increments  are  not  produced  for 
zero  deflections.  The  errors  introduced  by  this  assumption  are  currently  mitigated  by  the  robustness 
resulting  from  feedback  control  laws.  A  method  for  improving  the  performance  of  the  feedback 
control/control  allocation  system  is  presented  that  directly  attacks  the  inaccuracies  introduced  by  these 
linear  assumptions.  The  approach  is  demonstrated  using  a  dynamic  inversion-based  control  law  developed 
for  a  lifting  body  with  four  control  surfaces.  The  commanded  body  rates  form  the  input  to  a  dynamic 
inversion  control  law  that  forms  the  inner-loop  control  structure.  The  outputs  of  the  dynamic  inversion 
algorithm  serve  as  the  inputs  to  a  linear  programming  based  control  allocator.  The  control  allocation 
algorithm  is  a  mixed  optimization  scheme  that  minimizes  the  norm  of  the  error  and  the  deviation  of  the 
control  input  from  a  preferred  value.  Whether  the  underlying  dynamic  system  is  linear  or  nonlinear  in  the 
way  that  the  controls  affect  the  system  state,  the  control  allocator  requires  a  linear  approximation  to  the 
system  inputs.  In  this  work,  instead  of  using  a  linear  approximation  to  describe  the  control  variable  rate 
due  to  a  control  surface  deflection  (a  control  subspace,  that  contains  the  zero  element),  an  affine 
relationship  is  utilized.  This  allows  one  to  modify  the  input  to  the  control  allocator  by  providing  an 
additional  intercept  term  that  corrects  for  the  errors  introduced  by  the  control  allocation  algorithm’s 
assumption  of  linearity.  Since  on-line  control  allocators  perform  computations  at  flight  control  system 
update  rates,  the  decision  about  how  far  to  move  the  control  surface  in  the  next  time  interval  is  critical.  In 
order  to  accurately  compute  the  next  set  of  surface  deflection  commands,  the  local  slope  at  the  current 
operating  point  is  used  and  the  intercept  is  accounted  for  at  the  input  to  the  control  allocator. 

1.  Introduction 

The  forces  and  moments  produced  by  aerodynamic  control  surfaces  are  often  nonlinear  functions  of 
control  surface  deflection.  This  phenomenon  limits  the  accuracy  of  state-of-the-art  control  allocation 
algorithms  since  all  of  the  approaches  are  based  on  the  assumption  that  the  control  variable  rates  are  linear 
functions  of  the  surface  deflections  and  that  the  control  variable  rate  increments  are  not  produced  for  zero 
deflections.  The  errors  introduced  by  this  assumption  are  currently  mitigated  by  the  robustness  resulting 
from  feedback  control  laws.  A  method  for  improving  the  performance  of  the  feedback  control/control 
allocation  system  is  presented  that  directly  attacks  the  inaccuracies  introduced  by  the  linear  assumptions. 

Control  allocators  are  used  in  conjunction  with  some  type  of  feedback  control  law  whose  output 
consists  of  one  or  more  pseudo-control  commands.  The  number  of  pseudo-control  commands  is  always 
less  than  or  equal  to  the  number  of  control  effectors.  Dynamic  inversion  control  laws  [1]  and  control 
allocation  algorithms  fit  together  quite  naturally  since  the  pseudo-control  commands  are  easily  identifiable. 
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Dynamic  inversion  control  laws  require  the  use  of  a  control  mixer  or  control  effector  allocation  algorithm 
when  the  number  of  control  effectors  exceeds  the  number  of  controlled  variables,  or  when  actuator  rate  and 
position  limits  must  be  taken  into  account.  It  is  quite  common  that  the  desired  commands  can  be  achieved 
in  many  different  ways  and  so  control  allocation  algorithms  are  used  to  provide  unique  solutions  to  such 
problems.  Numerous  control  allocation  algorithms  exist  and  surveys  of  existing  techniques  are  presented  in 
Page,  et.  al.  [2]  and  Bodson  [3]. 

As  stated  previously,  state-of-the-art  control  allocation  algorithms  all  use  the  assumption  that  a  linear 
relationship  exists  between  the  pseudo-control-commands  (daes)  and  the  control  effector  displacements  ( 8 ), 
i.e., 


ddes=B6.  (1) 

In  the  one-dimensional  case,  this  is  equivalent  to  a  linear  subspace  containing  the  zero  element,  in  the 
6,daes  plane.  In  reality,  the  pseudo-controls  are  nonlinear  functions  of  the  surface  deflections  as  well  as 
other  parameters  P  that  typically  include  angle-of-attack,  sideslip,  and  Mach  number  and  can  be 
represented  by 

daes=g(P.8)-  <2) 

At  any  given  deflection,  one  could  find  a  B  for  which  B6  =  g(P,8) ;  however,  the  individual  control 

effectiveness  parameters  (slopes)  would  be  inaccurate. 

In  the  one-dimensional  case.  Equation  (1)  is  the  equation  of  a  line  passing  through  the  origin.  Figure  1 
shows  a  plot  of  a  typical  moment  generated  by  a  control  surface  as  a  function  of  control  surface  deflection, 
as  well  as  the  local  and  global  slopes  (bactuai  and  bi,  respectively)  at  the  current  operation  point  5i.  As 
shown  in  Figure  1,  for  the  one-dimensional  case,  Equation  (1)  utilizes  a  global  slope  approximation;  slope 
bi  is  used  as  an  entry  in  the  B  matrix,  while  the  true  slope  is  bactuai*  I^^  order  to  enforce  the  condition 
B8  =  g(P,8) ,  one  must  purposely  introduce  an  error  in  the  local  slope,  from  bactuai  to  bi.  This  error  is  what 

provides  the  zero  element  in  the  linear  subspace.  In  Figure  1,  P  in  G5(P,S)  is  set  to  Pq  to  indicate  that  a 
particular  operating  point  has  been  established. 


On-line  control  allocators  must  produce  control  effector  commands  at  the  flight  control  system  update 
rate;  therefore,  the  decision  about  how  far  to  move  the  control  surface  in  the  next  time  interval  is  critical. 
In  order  to  more  accurately  compute  the  set  of  surface  deflection  commands  at  the  next  time  instant,  the 
local  slope  at  the  current  operating  point  can  be  used  and  the  error  introduced  by  using  the  local  slope  in  the 
B  matrix  can  be  corrected  by  using  an  intercept  ( e(6) )  to  adjust  the  pseudo-control  commands.  In  other 

words,  we  make  use  of  an  approximation  for  g(P,8)  of  the  form: 
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(3) 


d<fcs=g(P.«)=B8  +  €(8). 

We  now  form  a  new  pseudo-control  command  d'des  such  that 

d'des=ddes-e(8)  =  B8.  (4) 

In  a  discrete-time  implementation,  one  would  use  the  following: 

^kdes  ^kdes  “  ® k^k+l 

where  k  indicates  the  k*  time  instant.  At  this  point,  any  of  the  current  linear  control  allocation  approaches 
can  be  used  to  approximate  the  exact  solution.  Equation  (2),  with  greater  accuracy.  In  this  work,  we  show 
that  utilizing  an  affine  relationship  between  the  pseudo-control-commands  and  the  control  effector 
displacements  can  yield  improvements  to  the  control  system’s  performance. 

Using  the  one-dimensional  case  as  an  example.  Equation  (3)  is  equivalent  to  an  affine  subspace,  that  is, 
one,  which  does  not  contain  the  zero  element.  Stated  in  another  fashion,  £(8)  represents  the  “y-intercept” 
term  of  a  line  in  a  plane.  For  visualization  purposes,  consider  a  vehicle  with  one  control  surface  and  one 
axis  to  control.  Figure  2  shows  the  effect  of  including  the  “intercept”  term  in  the  approximation  to  the 
control  inputs.  Contrary  to  the  situation  described  in  Figure  1,  in  this  case,  the  zero  element  is  not  required. 
Hence,  local  slopes  with  non-zero  G5(Po,8)  intercepts  are  acceptable.  The  parameters  bi  and  ba  represent 
the  slopes  of  the  lines  that  are  tangent  to  the  moment  curve  at  5i  and  82,  respectively.  In  other  words,  these 
slopes  are  the  local  control  derivatives. 
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Figure  3  shows  the  impact  of  e(5)  on  the  moments.  If  the  control  surface  is  operating  at  5i  without  Ei(5i), 
then  the  selected  moment  is  GsfPoiSi)  (see  dashed  lines).  However,  the  selected  moment  should  be 
G6(Po,5i)  +  ei(6i).  Hence,  with  the  additional  affine  term  included,  more  accurate  moments  are  generated. 

G  6(P  0,5) 

G6(Po,62)+  62(82 

G8(Po,5,)+  e,(6i 

G  8(P  0,8, 

G  8(P  0,82 


Figure  3.  Effect  of  e(8)  on  Moment  Selection. 

In  either  of  the  cases  discussed  above.  Equations  (I)  and  (3),  the  moments  generated  by  the  vehicle  are 
set  equal  to  either  a  purely  linear  term  or  an  affine  term.  Typically,  a  nonlinear  relationship  exists  between 
these  quantities;  however,  this  nonlinear  relationship  cannot  be  solved  with  any  control  allocation 
algorithm,  which  makes  the  assumption  that  the  controls  can  be  represented  by  a  linear  or  affine  term.  The 
best  one  can  do  is  to  reduce  the  errors  in  the  approximation,  hence  the  inclusion  of  £(8).  The  nonlinear 

behavior  is  a  direct  consequence  of  the  model,  in  particular,  how  the  inputs  affect  the  states. 

One  last  point  that  should  be  made  is  the  case  where  the  moment  curve  is  not  a  monotonically 
increasing  or  decreasing  function  of  5.  Assume,  for  example,  that  the  moment  curve  is  as  shown  in  Figure 
4.  This  control  surface  has  the  unusual  feature  that,  for  5  >  5p  (the  control  surface  location  where  the  peak 
moment  is  produced),  the  moment  generated  is  less  than  the  moment  produced  at  5  =  8p. 

G 

G  8(P  0.5 
G  8(P  0.5  „ 


T  ypical  M  om  ent 
Curve 

Figure  4.  Non-Monotonic  Moment  Curve. 

If  the  control  surface  is  operating  in  the  region  where  8  >  8p,  then,  when  a  smaller  moment  is  desired, 
instead  of  reducing  the  actuator’s  position  towards  zero,  the  control  allocation  algorithm  may  actually 
increase  the  actuator’s  position,  that  is,  farther  to  the  right  of  8p.  In  a  local  sense,  this  may  be  acceptable 
behavior,  but  due  to  actuator  position  limits,  there  is  a  point  at  which  8  =  8„a,.  The  control  allocation 
algorithm  will  work  correctly  as  long  as  the  desired  moment  is  between  G8(Po,8inax)  and  GglPcSp).  If  a 
moment  smaller  than  G5(Po,8max)  is  requested,  then  the  actuator  position  should  be  in  a  region  where  8  <  8p. 
However,  the  algorithm  may  not  move  over  the  peak  in  the  moment  curve,  effectively,  limiting  the 
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actuator’s  position  to  6p  <  5  <  8max.  This  is  similar  to  the  local  minimum  problems  that  occur  in  numerical 
methods.  Intuitively,  what  should  occur  in  this  situation  is  the  allocation  algorithm  should  move  over  the 
“hump”  in  the  moment  curve  to  the  region  where  5  <  6p.  Fortunately,  this  type  of  actuator  behavior  is  not 
extremely  prevalent  in  flight  control,  however,  if  it  does  occur,  care  must  be  taken  to  ensure  proper  results 
from  the  control  allocation  algorithm.  Two  solutions  to  this  problem  are  as  follows:  force  the  moment 
curve  to  a  monotonic  function  of  5  using  clipping  logic  (see  dotted  line  in  Figure  4)  or  reduce  the  position 
limits  on  5  to  a  point  where  the  adverse  sections  of  the  moment  curves  are  not  available  (5  <  5p).  One 
should  notice  that  altering  the  position  limit  of  the  control  surface  may  be  undesirable  because  the  adverse 
section  of  the  moment  curve  for  one  axis  may  be  beneficial  to  another  axis.  Likewise,  clipping  the  moment 
curve  results  in  inaccuracies  in  the  control  allocation.  Hence,  one  must  carefully  balance  the  trade-off 
between  the  two  options. 

Although  the  theory  described  in  this  section  was  demonstrated  using  one-dimensional  surfaces,  this 
approach  is  directly  applicable  to  higher-order  dimension  problems.  If  two-dimensional  (two  desired 
moments)  are  considered,  then  the  moment  curves  in  the  previous  figures  become  moment  surfaces,  which 
are  functions  of  two  independent  variables. 

2.  Model  Development  and  Dynamic  Inversion 

For  the  purpose  of  demonstration,  we  develop  a  dynamic  inversion  control  law  for  a  lifting-body 
reentry  vehicle  with  four  control  surfaces.  The  control  surfaces  include  two  flaperons  and  two 
ruddervators. 

An  outer-loop  control  system  generates  body-axes  angular  velocity  commands  (pdes,  qdes»  ^des)*  that  the 
inner-loop  attempts  to  track,  with  the  use  of  a  dynamic  inversion  control  law.  The  development  of  the 
vehicle  model  follows  the  work  of  Doman,  et.  al.  [4].  The  rotational  dynamics  of  a  lifting  body  can  be 
written  as: 

a)  =  f((0,P)  +  g(P,8)  (6) 

where  (o  =  [p  q  r]  and  P  denotes  measurable  or  estimable  quantities  that  influence  the  body-axes  states. 
The  parameter  P  includes  variables  such  as  angle  of  attack,  sideslip,  Mach  number,  and  vehicle  mass 
properties.  The  term  g(P,6),  in  Equation  (6),  includes  the  control  dependent  accelerations,  while  the  term 
f(a),P)  defines  accelerations  that  are  due  to  the  base-vehicle’s  aerodynamic  properties.  The  moment 
equations  for  a  lifting  body  in  body  axes  [5]  can  be  manipulated  to  form  the  control  dependent  and  control 
independent  terms  in  Equation  (6).  The  mass  properties  of  the  vehicle  under  consideration  are  constant, 

thus,  the  time  derivative  of  the  inertia  matrix  can  be  set  to:  1=0.  Equation  (6)  becomes: 


CO 


=  I  ^(Gb  -coxIco) 


where 


Gb=Gbae(®.P)+G8(P.8)  = 

l" 

M 

+ 

'l' 

M 

N 

BAE 

_N_ 

(7) 


(8) 


I  is  the  inertia  matrix,  and  L,  M,  and  N  are  the  rolling,  pitching,  and  yawing  moments.  In  Equation  (8), 
Gbae(®’P)  the  moment  generated  by  the  base  aerodynamic  system  and  G5(P,8)  is  the  sum  of  moments 
produced  by  the  control  effectors.  Therefore, 

f(0).P)  =  r‘[GBAE(“.P)-«>Xl®]  (9) 


5 


and 


g(P,8)=r*G8(P,6). 


(10) 


In  order  to  utilize  a  linear  control  allocator,  it  is  necessary  that  the  control  dependent  portion  of  the 
model  be  linear  in  the  controls.  Hence,  a  linear  approximation  is  developed  such  that: 

G8(P,8)  =  Gs(p)6+e(P,6).  (11) 

Using  eq.  (3),  the  model  used  for  the  design  of  the  dynamic  inversion  control  law  becomes 

m=f(a),P)+r‘G8(P)8+r's(P,8).  (12) 


The  objective  is  to  find  a  control  law,  which  provides  direct  control  over  co ,  so  that  ra  =  o>des  •  Hence,  the 
inverse  control  law  must  satisfy: 


Wdes-f(«>,P)-r‘8(P,8)  =  r‘G8(P)8 .  (13) 

Since  there  are  more  control  effectors  than  controlled  variables  and  the  control  effectors  are  restricted 
by  position  and  rate  limits,  a  control  allocation  algorithm  is  used.  For  the  lifting  body  under  consideration, 
there  are  three  controlled  variables,  namely,  roll,  pitch,  and  yaw  rates,  while  there  are  four  control  surfaces. 
Hence,  a  control  allocation  scheme  must  be  used  to  insure  that  Equation  (13)  is  satisfied.  The  control 
allocation  scheme  used  in  this  work  draws  heavily  on  the  work  of  Buffington,  et.  al.  [6,7]. 

There  are  two  branches  in  the  control  allocator,  a  sufficiency  branch  and  a  deficiency  branch.  The 
deficiency  branch  is  used  to  test  the  feasibility  of  Equation  (13),  while  the  sufficiency  branch  is  used  to 
optimize  a  sub-objective  when  excess  control  authority  is  available.  For  convenience,  let  Equation  (13)  be 
represented  by 


4s=^s-f(o),P)-8(P,8) =r'G8(P)S=B8-  (14) 

If  it  is  not  feasible  to  obtain  the  relationship  shown  in  Equation  (14)  due  to  control  effector  restraints,  then 
the  difference  between  the  desired  and  actual  effector-induced  body-axis  accelerations  is  minimized.  A 
complete  discussion  of  the  two-branch  control  allocation  scheme  can  be  found  in  Buffington  [6].  Equation 
(14)  is  in  a  form  that  makes  the  dynamic  inversion  more  accurate  since  local  slopes  can  be  used  to  populate 
the  B  matrix  and  the  errors  introduced  by  using  local  slopes  can  be  corrected  by  the  intercept  vector. 

3.  On-line  Computation  of  Slopes  and  Intercepts 

Aerodynamic  tables  that  can  be  used  to  compute  slopes  and  intercepts  for  the  problem  under 
consideration  are  generally  large  and  the  interpolation  process  can  be  slow.  A  common  solution  to  storing 
these  tables  is  to  fit  curves  to  the  tables  and  store  a  set  of  basis  functions  and  coefficients  so  that  table 
entries  can  be  estimated  on-line.  In  this  section  we  wish  to  point  out  that  the  slope  and  intercept  data  can 
be  stored  rather  compactly  using  least-squares  curve  fits  to  an  aerodynamic  table  that  contains  the  moment 
increment  due  to  an  effector  deflection  over  a  range  of  effector  independent  parameters  (P).  Let  g(P,8) 
denote  the  nonlinear  effector  data  stored  in  the  table  and  let  R  denote  a  set  of  basis  functions  that  can  be 
weighted  and  summed  such  that: 


g(P,6)  =  F'^R,  (5)  +  G'^R2  (P)  +  H'^R3(P,8)  (15) 

where  F,  G  and  H  are  vectors  of  constant  coefficients.  Thus,  a  curve  fit  for  the  elements  of  the  control 
effectiveness  matrix  (B)  can  be  obtained  from: 


6 


(16) 


a  g(P , «)  _  3  Ri(«) .  a  RjCP.  8) 

~Tr~ 

Therefore,  a  curve  fit  for  e(8)  is  given  by: 

e(^sF'^  R,(8)-^^1^8  +G’^R2(P)+H^  R3(P,6)-^^^^^---5  •  (17) 

d8  CrO 

Thus,  the  same  coefficients  that  are  used  to  store  the  slopes  are  also  used  to  store  the  intercepts. 

4.  Results 

A  six  degree-of-freedom  simulation  was  constructed  for  a  lifting  body  with  four  control  surfaces. 
An  inner-loop  control  system,  consisting  of  prefilters,  control  allocation  algorithm,  and  the  plant,  was 
constructed  for  the  purpose  of  tracking  commands  in  the  rolling,  pitching,  and  yawing  moments.  The 
inner-loop  control  structure  is  shown  in  Figure  5.  An  outer-loop  control  system  was  also  developed  to 
control  the  vehicle’s  flight-path  angle,  as  shown  in  Figure  6.  In  Figures  5,  the  notation  CV  is  used  to 
represent  command  variables. 


Figure  6.  Flight-Path  Angle  Controller. 
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Although  Figure  6  displays  the  proportional-derivative  flight-path  angle  controller  as  only  controlling  the 
pitch  axis,  some  other  minor  compensation  is  performed,  such  as,  regulating  sideslip. 

The  trajectory  that  was  used  in  the  simulations  is  as  follows:  the  vehicle  is  trimmed  at  an  initial  flight- 
path  angle  of  -27".  At  time  =  2  seconds,  a  -50"  flight-path  angle  is  commanded.  ^When  time  =  13  seconds, 
the  vehicle  pulls  up  to  -27"  then,  at  time  =  23  seconds,  the  vehicle  flares  to  a  -2.5"  flight-path  angle.  Figure 
7  shows  a  time  history  of  the  actual  glide  slopes,  with  and  without  £(8)  compensation,  and  the  commanded 
glide  slope.  Here,  it  can  be  determined  that  the  actual  glide  slope  tracks  the  command  more  accurately 
when  £(8)  is  used. 


Figure  7.  Flight-Path  Angle  Command  and  Responses. 

In  order  to  examine  the  effectiveness  of  the  additional  control  allocation  compensation,  it  is  pertinent 
to  consider  the  desired  pitch  rate  and  the  actual  pitch  rate.  Figure  8  shows  the  error  in  the  pitch  rate  from 
the  glide  slope  control  simulation  run.  The  two  curves  in  Figure  8  represent  the  errors  in  pitch  rate  with 
and  without  e(8)  compensation.  The  solid  line  represents  the  pitch  rate  error  with  8(8),  while  the  dotted 
line  shows  the  pitch  rate  error  without  8(5) .  It  can  be  seen  that  the  pitch  rate  errors  are  significantly 
reduced  when  8(6)  is  included  in  the  control  system. 


Figure  8.  Error  in  Pitch  Rates. 
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4.  Conclusions 


In  a  linear  programming  control  allocation  algorithm,  the  systems’  typically  nonlinear  input  is  forced 
to  be  equal  to  a  linear  expression.  Hence,  fitting  errors  certainly  affect  the  preciseness  of  the  control 
system.  To  help  improve  the  fit,  an  additional  term  was  identified  and  added  to  a  control  allocation 
scheme.  This  term  has  the  effect  of  allowing  the  desired  set  of  body-axis  moments  to  be  represented  by  an 
affine  subspace,  instead  of  a  linear  subspace.  In  this  way,  two  terms  are  required  to  represent  the  desired 
set  of  moments,  a  local  slope  and  a  corresponding  “intercept”  term.  Some  of  the  benefits  of  including  this 
term  in  inner-loop  control  systems  were  also  highlighted. 
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